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ABSTRACT 


Methods are presented for calculating design limit loads compatible with 
probabilistic structural design criteria. The approach is based on the 
concept that the desired "limit load," defined as the largest load occurring 
in a mission, is a random variable having a specific probability distribution 
which may be determined from extreme-value theory. The "design limit load," 
defined as a particular value of this random limit load, is the value 
conventionally used in structural design. Methods are presented for deter- 
mining the limit load probability distributions from both time-domain and 
frequency-domain dynamic load simulations. Numeriial demonstrations of the 
methods are also presented. 
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l.O INTRODUCTiON 

The purpose of this report is to describe and numerically demonstrate methods 
for combining payload parameter variations with the input environment in 
probabilistic structural design loads analyses. The design loads resulting 
from these methods are compatible with probabilistic structural design criteria. 
The approach is based on the concept that the desired "limit load," defined 
as the largest load occurring in a mission, is a random variable having a 
specific probability distribution v/hich may be determined from the extreme- 
value theory of probability. The "design limit load," defined as a particular 
value of this random limit load, is the value conventionally used in structural 
design. 

The scope of this study was limited in three general areas. First, no attempt 
was made to include the effects of structural fatigue. The technical theory 
is concerned only with structural designs corresponding to the single applica- 
tion of an extreme load to an undamaged structure. Second, no attempt was 
made to define rationale for selecting acceptable probabilities of failure 
to be used in the structural design criteria. Third, the technical theory 
is concerned only with the preliminary design/redes ign/des ign verification 
phases of a project. No attempt v<as made to address the inverse problem of 
operational constraints and decisions. 

A discussion of a proven general probabilistic structural design approach 
is presented in Section 2.0 along with some, basic results of extreme-value 
theory which are particularly applicable to structural loads. Section 3-0 
presents methods for determining extreme-value limit-load probability 
distributions from conventional time-domain and frequency-domain dynamic 
loads analyses. Numerical demonstrations of each of these methods are 
presented in Section k.D. Conclusions from the present research and 
recommended a.'eas for future research are presented in Section 5*0. A 
comprehensive list of references completes this report. 
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2.0 THEORETICAL BACKGROUND i 

1 

The concept of a randomly varying limit load described by a theoretically 
correct p.-obability distribution and the use of a particular value of 
this random limit load for structural design purposes are of basic import- 
ance in probabi 1 istic structural design criteria. Since the limit load 
is conventionally defined as the largest load occurring in a mission, 
the probability theory of extreme values is useful in determining the 
theoretically correct limit-load probability distribution. Section 2.1 
contains some basic results of extreme-value theory which are particularly 
applicable to structural loads. Since the determination of probabilistic 
structural loads is meaningful only within the larger context of structural 
design. Section 2,2 includes details of the application of probabilistic 
load quantities in a general structural design approach. 






2.1 Limit-Load Probability Distributions I 

The limit load for a structural component is conventionally defined as the 
largest load occurring durini, =» given mission. The probability that the 
component load x is the largest value among n independent observations is 
defined by 

^„(x) = If(x)3" (1) 

where F(x) is the underlying cumulative distribution function (CDF) for the 
load. Thus is, by definition, the cumulative distribution function 

of the limit load for a mission v^hich has n Independent occurrences of applied 
load. The probability theory of extreme values, as presented by Gumbel 
(Reference 1), is concerned with describing the limit-load distribution 
function (l>^) for various forms of the underlying distribution (F) . 


Two parameters frequently used in extreme-value theory are the characteristic 
largest value and the extremal intensity function. The characteristic 
largest value (u^) in a Sample of n observations is defined by Gumbel 
(Reference 1, page 82) in terms of the following equation: 


v/here is the underlying CDF evaluated at the characteristic largest value 

Thus, as indicated by Equation (2), u^ is that value of the random variable 

v;hich v/ill be equalled or exceeded one time in n observations, on the average. 

The extremal intensity function (a ) in a samole of n observations is defined 

n 

by Gumbel (Reference 1, page 81») as follows: 


f(Un) 

“n = T-"f(u')‘ 
n 


(3) 


where is the underlying probability density function (PDF) and 

the underlying CDF, both evaluated at the characteristic largest value. The 
inverse of the extremal intensity function, called Mill's ratio, is tabulated 
by K. Pearson for the normal distribution (Reference 2, page 11). 
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The underlying distribution F(x) is said to be of the exponential type if 
f(x) approaches zero for large [x| at least as fast as the exponential 
distribution, f(x) = \e . For any distribution of the exponential type, 
Gumbel (Reference 1, page 16F) shoi-/s that the CDF for large x is aoui'ox- 
imately equal to 


-a {x-u ) 
F(x) = 1 - S. " ” 




An asymptotic distribution of extrene largest values can be obtained by 
substituting Equation (A) into Equation (1) and taking the limit as n 
becomes infinite 


4-^^’^{x) = £ 

n-w 


-a (x-u ) 
. , e n n 

im I - — 
n 

L 


(5) 


Evaluating this limit by means of the logarithmic series results in the 
first asymptotic distribution of extreme largest values, subsequently 
called the extremal type 1 distribution: 


^ (x) = exp 




The corresponding PDF, v^hich is positively skewed, is given by 


. 0 ) 


-a (x-u ) I 

(x) = a exp -a (x-u ) -e j 

n ^ n n | 


( 6 ) 


(7) 


The most probable value or mode (m^) of this distribution is equal to the 
characteristic largest value; 


m = u 
o n 


( 8 ) 


The fifty-percentile value or median (n ) is given by 


ini-ZnO-S) = u + 0.36651292 

m = u - n 

e n a a 


(9) 


The mean (m) is 


m = u + E 
n — 
a 


where = 0.5772156b is Euler's constant. 


Oo) 



The standard deviation (s) is given by 


s = 


/k a 


and the coefficient of variation (V = s/m) is 


V = 


/6(a u + C,.) 
n n E 


i 


I 


( 11 ) 


( 12 ) 


Equations (10) through (12) define parametric values for the extremal tyoe 
I distribution corresoond ing to a single mission. Parametric values for 
the largest load occurring in N missions are as follows; 


^ 0 « 
m., = m + — JlnN* 

N 7T 


(13) 


s 


= s 


(Ilf) 


" (1+1^ inU-y) 

Tf 

These relations are derived in Reference 3 (page 67)- Note that the 
standard deviation (s) and the extrema! int-.'.nsity function (:i^) for the 
extremal type I distribution are theoretically independent of sample size. 


According to Gumbel (Reference 1, page 182) the extremal type I distribution 
is often satisfactorily represented by the lognormal distribution. The 
lognormal distribution with coefficient of variation equal to 0.3olf is 
essentially Identical to the extremal type I distribution. For coefficients 
of variation betv/een 0.31 and O.UZ, the extremal and lognormal distributions 
are graphically indistinguishable. An example of the validity o'' the log- 
normal appro-ximat ion to the extremal type I distribution is given in Reference 
A. For this analysis, 28 sets of internal lead quantities v/ere calculated 
as the maximum values experienced in each of 100 simulated lunar landings. 

A Chi-square test of the hypothesis that the loads were lognormal ly dis- 
tributed resulted in cumulative probabilities ranging from 5 to SO percent. 

The iognormal approximation was therefore considered acceptable since the 
Chi-square probabilities v/ere less than 90 percent for all 28 internal load 
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quantities. The Chi~square hypothesis is usually accepted for cunulative 
probabilities as high as 99 percent. The coefficients of variation for 
these load quantities varied between 0.2 and 0.4. 

The extrenal type I distribution, defined by Equations (6) through (15), 

Is the theoretically proper distribution for limit loads due to any con- 
dition having an exponential-type underlying probability distribution and 
a sufficiently large number of independent load occurrences. For the 
exponential distribution, convergence to the asymototic extremal tyoe I 
distribution Is essentially complete for 100 observations (Reference 1, 
page 116). For the normal distribution, however, convergence to the 
asymptotic type 1 distribution is extremely slow. According to Fisher and 
Tippett (Reference 5, page 189), close convergence is attained only for 
sample sizes on the order of 10^^- Such large samples correspond to 
characteristic largest values of the standardized normal variate on the 
order of 16. 


Accurately describing extreme values from an underlying normal distribution 
is necessary due to the central role of the normal distribution in engineer- 
ing applications. The theoretical distribution of extreme largest values 
from variously sized samples of standardized normal variates was tabulated 
by K. Pearson in Reverence 2 (page 162). Plots of these tabulated values 
on lognormal probability paper indicate that, for certain sample sizes, the 
theoretical distribution of normal extremes can be adequately approximated 
by cne lognormal probability distribution. In fact, the theoretical dis- 
tribution plots essentially as a straight line on lognormal probability 
paper for standardized characteristic largest values (Uj^) of approximately 

2.16. This value of u corresponds to a sample size (n) of approximately 65. 

n 

The lognormal PDF may be written for the normal extreme variate (x) as 
fol lov/s : 


f (x) = exp 

/ItT 6x 


1 f ilnx- y 

2 ^ 6 ' 


(16) 


v/here y 
6 


is the mean of £nx and 

is the standard deviation of Hnx. 
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The parameters y afid 5 used for the lognormal approximation are obtained 
in terms of the standardized extreme median (x) and standardized extreme 
mode (k) by means of the follov/ing identities: 

(1) The mean of £nx is the logarithm of the median of x; 

( 2 ) The variance of £nx is the logarithm of the ratio of 
median of x to the mode of x (Reference 1, page 18). 

Let the underlying normal distribution of interest have mean U and standard 
deviation g. The required lognormal parameters are then given by 


Y = i.n(p + xo) 



07 ) 

08 ) 


The median (x) of the standardized normal extreme for n samples is defined 
b'' the follov;ing equation; 


IfOOj" = 0.5 09 ) 

Combining Equations (2) and (1?' to eliminate n gives the following desired 

*J 

equation for the standardized extreme median (x) in terms of the standardized 
characteristic largest value (u^) : 

F(x) = exp [(in 0 . 5 ) (I " F(u^))2 (20) 

where F is the normal CDF. 


The mode (x) of the stan. ardized normal extreme for n samples is defined by 
Gumbel (Reference 1, pagi 133) in terms of the following equation. 

x«F(x) = n-1 ^ 2 J) 

f(x) 

where F is the normal CDF, and 
f is the normal PDF. 

Combining Equations (2) and (2l) to eliminate n gives the following desired 

equation for the standardized extreme mode (x) in terms of the standardized 

characteristic largest value (u ): 

n 
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x.F(x) 

f(x) 


l-F(u_) 


( 22 ) 


The lognormal approximation to the distribution of normal extremes is 
defined by Equations (16) through (22). ~his representation may be 
considered adequate for values of the standardized characteristic largest 


value (u ) less than 3- 
n 


A second approximation to the distribution of normal extremes was proposed 
by Fisher and Tippett (Reference it). The proposed CDF is of the form 


$^3)(y) ^ 


r 

-(■^) 

u 

^ n -* 


where k » 


(uj - 1) 


(23) 


(2b) 


Tnis general form is denoted by Gumbel (Reference 1, page 258) as the 
third asymptotic distribution of extreme values or the extremal type III 
distribution. By inverting Equation (23), approxi-mate percentage points 
for extremes of the standardized normal variate are obtained as follows 
in terms of the cumulative probability, p: 

X - u _ Ihl'^np) 


exp 


[ , ' in(-£np) 

,.n u ; 

n k 

J 


(25) 


A special characteristic of this extremal type 111 distribution is that it 
converges for increasing values of the parameter k toward the extremal type 
I distribution. Thus, in practice, the extremal type III distribution may 
be used to represent normal extremes for all values of u^ greater than 3- 

For very large values of u^, the extremal type I distribution, which is the 
theoretical asymptotic distribution of normal extremes, may be used- For the 
normal standardized variate (y) , the type 1 distribution function is 

4 ,”'(y) = exp(-e 




(26) 


where 


.. . . . -c 

a = u + u -2u + lOu 

n n n n n 


(27) 
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This expression for the standardized normal extremal intensity function 
derived by Gumbel (Reference 1, page 137)* The type I extreme value dis 
tribution may be used, if desired, to describe normal extremes for 
standardized characteristic largest values exceeding 8. 
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2.2 


Probabilistic Structural Design Criteria 


i 

The extended reliability structural design approach proposed by Ang and 
Atnin (Reference 6) recognizes both the probabi 1 istic nature of limit loads 
and strength and the analytical uncertainties associated with their evalua- 
tion. The uncertainties associated vjith determining limit loads and strengths 
can be quantified by a factor of uncertainty (v) equal to or greater than 
unity. Thus the event (S/L ^ v) constitutes a state of structural safety, 
where S and L represent the strength and limit load associated with a struc- 
tural component. If S and L are both random variables, then the probability 
Pis /L ^ v] is -p proper measure of structural safety. The extended reliabil- 
it> structural design approach is then expressed by the following probabil- 
isllc equation for structural safety: 

Pi S ^ v] = } - Pp (28) 


where S is tha random variable describing the component strength, 

L is the random variable describing the component limit load, 

V is the Ang-Amin coefficient of uncertainty for the component, and 
Pp is the component probability of failure or acceptable risk. 


When the limit-load and strength probability density functions are known, 
the structural design approach may be expressed in tv;o equivalent forms; 


1 


y/v 


f|^(x) • fj(y) dx dy 


(29) 




vx 


fs(y) • 


(30) 


where fj^(x) is the limit-load probability density function (PDF), and 
fg(y) is the strength PDF. 


The conventional factor of safety is defined as the ratio of the allowable 





FS - ^ (31) 

" 1 

where is the value of the random strength corresponding co a 
specified exceedance orobability i 
Ljj is the value of the random limit load corresponding to 
a spec if iced non-exceedance probability (Pp) • 

The purpose of the factor of 'afety in the structural design procedure is 
to locate the strength POP alative to the given limit-load PDF so that 
Equation (29) or (30) results in the required component probability of 
failure. This concept is illustrated in Figure 2.2-1. For most probability 
distributions, the integral of Equation (29) or (30) must be evaluated 
numerically and the required factor of safety determined by tr ial -and-error 
procedures. However, for certain specific distributions, closed-form evalua- 
tions leading to convenient design formulas are possible. 

A particularly convenient design factor-of-safety equation occurs when both 
limit loads and strengths are assumed to follov; the lognormal probability law. 
As discussed in Section 2.1, the lognormal distribution often accurately 
represents the theoretically proper distribution for limit loads. Moreover, 
for much existing strength data, the lognormal distribution also is a 
satisfactory representation, due perhaps to the deletion of low-strength 
values by quality-control procedures. 



The component factoi — of-safety expression for lognormal limit loads and 
strengths is derived in Reference 3 in the follov;ing form; 

FS = 

+ F 

where v is the Ang-Amin coefficient of uncertainty, 

Pp is the probability of failure or acceptable risk, 

Pp is the non-exceedance probability for design limit load (S-jj) , 
P^ is the exceedance probability for allov/able stress (S^) * 
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PROBABILITY DENSITY FUNCTIONS 




Figure 2.2- 1. Graphical Representation of Factor of Safety 





and Vj are limit-load and strength coefficients of variation, 

I 

P (P) is the inverse of the standardized normal cumulative 


distribution function gi\en by 



rP"’(P) 

1 

/27 


it^ 

2 ^ dt 


(33) 


The numerical behavior of the lognomal/lognornal factor of safety is shown 
graphically in Figure 2.2-2. For this plot, the defining probabilities for 
design limit load and allowable strength are both taken as 99 percent, and 
the coefficient of uncertainty is taken as unity. The factor of safety is 
seen to increase monotonically with decreasing probability of failure for 
given load and strength coefficients of variation. 


From Equation (32), the component factor of safety corresponding to a 
specified probability of failure may be computed. The allowable strength 
is then determined, from Equation (31), as the product of the factor of 
safety times the design limit load. Additional details regarding the 
application of this probabilistic design approach are presented in Ref- 
erence 3* Procedures for determining the basic limit-load probability 
distributions from which the specific design limit load is selected are 
discussed in the following section. 
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3.0 


METHODOLOGY DEVELOPMENT 


Structural dynamic analyses resulting in design limit loads may be performed 
either in the time domain or in the frequency domain. The Taylor's series 
method and the Monte Carlo method are tvi/o widely used techniques for determining 
limit loads from time-domain analyses. The Taylor's series method, described 
in Section 3.1, is an extension of the parameter variation study often performed 
to evaluate sensitivity to parametric data uncertainties. The Monte Carlo 
method, described in Section 3-2, is a simulation of the loading condition 
using a random combination of vehicle parameters and environments. For each 
load quantity of interest, the maximum vaije occurring in each simulated 
mission is identified and recorded. The maximum load data from a number of 
simulated missions approximates the desired extreme-value limit-load 
distribution. 

In Section 3.3« a new method is presented for determining the extreme-value 
limit-load distribution from a frequency-domain analysis. This method 
determines the probability distribution of the extreme largest load value, 
for a stationary Gaussian random process, occurring within a given mission 




3-1 Limit Loads from Taylor's Series Analyses 


;■ A detailed discussion of the use of the Taylor's series method to estimate 

; limit-ioad probability distributions for aerospace launch vehicles is 

presented by Lovingood (Reference 7)- This application involves first 
analytically simulating the structural loads and responses encountered by a 

nominal launch vehicle flying through a moderately severe synthetic wind 

profile. The resulting loads are considered to be the nominal or mean 
values for the limit loan probability distribution. The peak or design 
limit load values, which are defined as the "3*o" values having non-exceedance 
' probabilities of 0.9987, are next obtained by computing the variations in 

load due to 3"0 variations in the significant vehicle parameters, taking 
\ the root-sum-square variations of each load quantity, and adding these to 

the corresponding mean values. 

\ This method is useful for efficiently predicting preliminary and interim 

\ 

Structural design loads. However it has the disadvantage of requiring a 
synthetic wind profile defined such that the mean values of all the limit 
, loads of interest are produced by the analytical simulations. Besides the 

difficulty of defining this proper synthetic environment, the Taylor's 
series method is based on three fundamental assumptions v;hicti may not be 
valid for particular appl ications. These assumptions will be discussed in 
fie brief derivation which follows. A similar derivation in Reference 7 is 
somewhat more detailed. 

The distribution of a nonlinear function of several random variables may be 
obtained by approximating the desired function as a linear function in the 
region of interest. The mean and standard deviation of a linear function 
of several independent random variables are knov-/n from elementary probability 

theory (Reference 8, page 1*8). If X., X , .... X are independent random 

' ^ ^22 2 
variables having means mj . m^ and variances Sj f s^ ^n ’ 

respectively, and if aj, a^, are constants, then a linear random 

function may be defined as follov/s: 
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- a,X, f ax ^ ( 3 « 

The mean of f is 

•= Q-^TA-^ ^ J'/?„ (35) 

Thus the mean of a linear combination of random variables is equal to the 
linear combination of the means. This result is valid even if the X's are 
dependent. 


The variance of f is 

' -f-a^ S,, 06) 

Thus the variance of a linear combination of independent random variables 

is equal to the sum of the products of variances and squared constants. In 

addition, if the X's are normally distributed, then f is also normally 

2 

distributed v/i th mean and variance s^; . 

A nonl inean function may be expanded in a Taylor's series about any given 
point as follows (Reference 8, page 62): 


/a,, /x, X„) = "V 




UK - ■ 

(A/7' hy 

V AM / 


+ hiOftfr or., er {-iy) 




If the higher order terms are negligible, the mean of f is, according 
to Equation (35), approximately equal to: 
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( 38 ) 








^ sS ^ 



are independent, the variance 

imately 

equal to 

-)f\ 




'E 

-1 


, ”‘j 


ing 


(39) 


Futhermore, if the X's are normally distributed, f is approximately normally 
distributed. If the X's are normally distributed and if the function is 
linear so that Equation. (37) contains no higher order terms, then the mean 
and variance are exactly as given by Equations (38) and (39) and the 
theoretical distribution of the function is the normal distribution (Refer- 
ence 9. page 90 ). 


The three assumptions in the use of the Taylor's series method are as 
follows: 

(1) that the higher-order terns in the Taylor expansion are negligible 
compared with the first-order terms, 

(2) that the X's are independent, and 

(3) that the X's are normally jistributed. 


The accuracy of design limit loads determined by the Taylor's series method 
depends in part upon how well the pa’^ticular physical simulation is 
represented by these three assumptions. In Section i*.l, a discussion of 
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the effects of these assumptions is presented along with numerical 
demonstrations of the method. So long as the potential disadvantages 
of this method are recognized, it remains an efficient and useful tool 
for estimating preliminary and interim design limit loads. 
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3-2 Limit Loads from Monte Carlo Simulations 

I The Monte Carlo method is a powerful and general tool for predicting 

' structural design loads. The method has been gaining wider acceptance for 

dynamic load studies of aerospace vehicles (References 4, 10, 11. 12). For 
this application, the method consists essentially of simulating a random 
loading phenomenon by combining deterministic and probabilistic variables. 

The limit-load probability distribution for each load quantity is then the 
distribution of the largest loads occurring in each simulated mission. For 
the launch vehicle load simulations described in Reference 10, the determin- 
istic variables included such vehicle parameters as mass and geometry, 
structural dynamic characteristics, propellant slosh parameters, and control- 
system parameters. The probabilistic variables for this study were restricted 
to descriptions of the wind environment. The wind was represented both by 
detailed measured wind profiles including turbulence and by filtered measured 
j wind profiles with the turbulence considered separately using power spectral 

density (PSD) methods. 

In general, probabilistic variables may include any factors not determin- 
istically known, including initial conditions, propulsion characteristics, 
alignment tolerances, and mass properties. For time-domain simulations, 
sample values of individual randam variables r>iay be generated using digital 
random number generators such as those described in References 13 and l^». 
Sample time histories of random processes such as wind turbulence can be 
generated from PSD data using the technique described in Reference IS- Of 
course, actual sample values or sample time histories from test data may be 
used directly as the random inputs to a Monte Carlo time-domain simulation. 

A major consideration in the general application of the Monte Carlo method 
is to reduce the required cost of simulation as much as possible. In 
Reference 16 (page 1^6), H. Kahn describes several such techniques. Two 
of these (Russian Roulette and Use of Expected Values) have been used 
I successfully in structural load analyses. Russian Roulette involves 
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concentrating the computational effort on cases of special interest. For a 
landing dynamics analysis, the cases of interest may be those having the 
largest initial kinetic energy which therefore result in the largest structural 
loads. For a flight loads analysis, the cases of interest nay be those 
having the wind profiles resulting in largest loads; the critical profiles 
are identified using very greatly simplified flight simulations. These cases 
identified as being of special interest are then analyzed using the more 
detailed simulation methods. The Use of Expected Values is merely a separation 
of computational tasks into what can be efficiently calculated analytically 
and what must be simulated by Monte Carlo methods. An example of this 
technique is the separation of the wind profile into small-scale turbulence 
(efficiently treated by PSD methods) and large-scale variations as described 
in Reference 10. 

Another technique which may be used successfully for determining probabilistic 
design limit loads is the statistical estimation method. As an extension of 
the normal confidence limit concept, this method is based on the generally 
valid representation of random limit loads by the lognormal probability law. 

The expression of the one-sided normal confidence limit as derived in Reference 
is valid for Monte Carlo samples of 50 or more observations. This expression 
can be simply modified as follows to be valid for samples as small as 20 
observations. 

Let Yj , ^2 Yp independent observations of a normal random variable 

with mean and standard deviation s^. The unbiased estimates of the sample 
mean and variance, v;hich are stochastically independent, are given by 
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According to Wilks (Reference 17. page 208) , the sample mean (m*) is normally 
distributed with mean (ra^) and standard deviation (s^/n) and the sample 
variance is distributed as follows 




(^»2) 


(rt -0 ~ 

2 . 

The Chi-square distribution with k degrees of freedom (xj^) is approximately 
normal for large k (Reference 17. page 183). However, a much more rapidly 
converging approximation is given by Boivker and Lieberman (Reference 8, 
page 556) : 


zx* A/ (\Z!c-i , l) 


m 


The close convergence of this approximation for 20 degrees of freedom is 
shown in Figure 3-2-1. Combining Equations (^2) and (^3) results in the 
following approximate distribution for the sample standard deviation for n 
as small as 20: 






Define the true a x 100 percentile load by 
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The statistical estimate of F is 
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From Equations (35), (36), and (Ai<) , the mean and variance of F are 


e[F] 3 ,,, 

The one-sided confidence limit equation is 
Equation (49) implies that 

- - fv< ■ £ f p] 


(47) 

(48) 


(49) 


(50) 


H Var[Fj 


where 


^ - 
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^ Jt 


-K(2 


Substituting Equations (45), (47), and (48) into Equation (50) and solving 
for the appropriate root of K yields 


K - S>-^A/n 


(51) 




where 


5 

C 


C - K^/tLn-l) 


Equation (51) niay be used with Equation (46) to determine the one-sided 
confidence limit for any probability level (a) and confidence level (6) 
so long as the sample size (n) is at least 20. 
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The statistical estimation method for jse i'l estimating the limit load 
probability distribution from at least 20 unbiased Monte Carlo observations 
Is described as follows for each load qujntity of interest: 

(1) Calculate the sample mean and standard deviation of the natural 
logarithms of the observed loads using Equations (1»0) ar •' (4l). 


(2) Calculate the one-sided confidence limits for several different 
probability levels (a) for a given confidence level (S) using Equations 
{kG) and (51) 


( 2 ) 

F Uzj /J 


m/ ^ KUij) ■ s* 

my + !<■ Ui, iS) ■ 

• • • 


(3) Solve for the mean and standard deviation of the logarithms which provide 
the least-squares fit to the follovjing equations: 

= %(/?) + Kui 

F {(<■'>■ ■ >riy(p + 


w Convert my(2) and 5^(6) to lognormal mean '"^(6} and coefficient of 
variation using the following standard expressions ^ 

Vy(^) - [ ex/? (Sy(?)^) - ij ( 52 ) 

Wx C( 2 ) - [ I ''x (>^y t«) 

Equations (52) and (53) are consistent with the following notation: 

Y - ^ X 

where X ^ /? O f q/ C Va* “ Sk X ) 

Y ^ NhrrAol ( my^ ^/) 


25 




1 

t 

This statistics] estimation method provides conservative estimates of the 
lognormal parameters of limit loads determined from at least 20 Monte Carlo 
simulations. The degree of conservatism in the estimated parameters is, 
of course, dependent on the confidence level (£) chosen. The estimated 
parametric values are also somewhat dependent on the particular probability 
levels (ai) chosen for the least-squares fit. A numerical demonstration 
of this method is presented in Section k.Z. 
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3-3 Limit Loads from Frequency - Domain Simulation 

t 

In the practical solution of random vibration problems, the dynamic char- 
acteristics of a structural system are usually assumed to be linear and 
deterministic, and the excitation is assumed to be random. Furthermore, 
the random excitation is usually assumed to be stationary and Gaussian 
with zero mean value, since the random process for the response can then 
be completely characterized by its power spectral density function (Reference 
page 89 ). Solutions to two random vibration problems for this special 
case of stationary Gaussian response are available in the literature (Ref- 
erence 18, page 293). The threshold-crossing problem is concerned with the 
expected rate at which a random process X(t) exceeds a certain value. The 
peak-distribution problem is concerned both with the probability distribu- 
tion of peak magnitudes in X(t) and with the expected rate of occurrence of 
the peaks. However, neither of these available solutions provides the 
extreme-value probability distribution required for probabilistic ultimate 
strength design. The objective of the present study is to determine the 
probability distribution of the extreme largest value, for a stationary 
Gaussian random process X(t), occurring within a given mission length. This 
required limit-load probability distribution will be expressed in terms of 
the power spectral dens i ty f uncti on (PSD) of the calculated load. 


The real autocorrelation function associated with a real-valued stationary 
random process X(t) may be defined by 


„ / , _ Jim 1 

^xx^"'^ ■ T-x» 2T 


X(t)*X(t+T) dt 


'-T 


(5*<) 


Equations relating the autocorrelation function and the power spectral 
density function (PSD) are known as the Wienei Khintchine relations (Ref- 
erence 19 , page 579). For a real-valued random process, such as the random 
load in a structural member, the defining equations may be written 
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as follov/s: 
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RCt) 


.09 

G (u) COSWT dw 

0 


(55) 


G(u) = — 


r°> 


R(t) coscjT dt 


(56) 


where G(;o) is the load PSD with frequency (^) in radians/second. 


The load PSD njay alternatively be written with frequency in Hz as follov/s 


R(i) 


rCf) cos2'.rfT df 

■ 0 


r(f) = k 


R(t) cos2-fT dx 


(57) 

(58) 


where f = u/2- in Hz. 

For some applications, the load PSD may be more Con>/eniently defined in 
terms of spatial frequency (radians per unit distance) and spatial distance 
instead of circular frequency (radians per second) and time. Equations (55) 
and ( 56 ) with appropriate notation changes may be used as the defining 
Wiener-Khintchine relations for such applications. 


With no loss or generality, a stationary random process may be assigned a 
zero mean value. The variance of such a real-valued random process is 
obtained from Equations ( 55 ) and (57) by evaluating the autocorrelation 
function for zero time lag. 


R(0) = 


G(o))d.i = 


r(f) df 


(59) 


Equations (5^) through (59) form a consistent set of definitions for use in 
harmonic analysis of stationary random processes. Since many authors use 
alternate forms of the Wiener-Khintchine relations (Reference 19 , page 5S0) , 
special care is required v/hen applying formulas for random vibration analysis. 
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Standard methods are available for computing the PSD of loads in a linear 
structure due to stationary Gaussian excitation (References 18, 20, 21). 

The output response PSD for the r'" calculated load quantity is given by 
the following general equation: 

G^Cw) = L>-ri*0uj)J lGj:(ja))]{L.^(ju)} (60) 

where {L. (ju)} is the column matrix of complex frequency responses 
^ ^ th 

for the r load quantity and for i excitation points, 

[jL^j*(ja)J is a row matrix of the complex conjugates of L.^(ju), 
and 

iG^(ju)] is the PSD matrix of input power spectral density 

functions for each of the i excitation points and cross- 
power spectral densities between the excitation points. 

The following development converts the Gaussian load PSD typically defined by 
Equation (60) into an extreme-value limit-load probabi 1 i ty d istr ibut ion 
required for probabilistic structural design. 


The critical parameter in the three distributions used for describing 
extreme normal variates is the characteristic largest value (u) . Its 
magnitude increases v/ith sample size until, as n becomes very large, it con- 
verges to the most probable value (mode) of the asymptotic extremal type I 
distribution (Reference 1, page 172). However, as described in Section 2.1, 
the convergence of the normal extremes to the type I distribution is so slow 
that the lognormal and extremal type III distributions must be used for small 
and moderately sized samples. The following development is based on expres- 
sing the characteristic largest value in terms of Rice's theorem for the 
expected number of threshold crossings per unit time. 


According to Rice (Reference 22, page 192), the expected rate of zero 
crossings from below for a stationary Gaussian process with zero mean is 
g i ven by 


E(N^(0) = 


■'0 


r r(f) df 




r(f) df 


(61) 
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With the PSD defined in radians per second (or radians per unit distance). 
Equation ( 61 ) becomes 

ElN^(O)] = ^ 

With the random structural load PSD defined by Equation (60) the integrals 
of Equations (61) (62) will converge vihenever the input PSD has a finite 

variance. 

The equation for the expected number of times per unit time or distance that 
the Gaussian load passes through the threshold value (0 with positive slope 
is given by Rice (Reference 22, page 192) as follows: 

ElN. (fJ3 = E[N (0)3 exp ( -\ ) ( 63 ) 

2o 

where E[N (0)3 is defined by Equation (61) or (62) and 

2 

q is defined by Equation (59). 

Equation ( 63 ) may also be found in Reference )8(page 297), Reference 2 ] 
k2) , and Reference 23 (page 5-121) among many other sources. It is restricted 
to stationary Gaussian random processes having zero mean values. Since the 
Gaussian model is commonly used to represent inflight atmospheric turbulence 
(Reference 23 , page 5-H6) and transonic buffeting (Reference 20), this 
restriction is not significant to most current engineering applications. 

The expected number of threshold crossings in a given time or difiance 
interval (T) is obtained simply by modifying Equation ( 63 ) as follows: 

,2 

EIN^(c)3t = T-EiN (0)3 exp( ) ((>^) 

* ' 2o 

where T defines the length of a mission. 

The desired characteristic largest value in a sample of size n, u , is 
defined as follows by Gunbel (Reference 1, page 82): "In n observations, 

the expected number of values equal to or larger than u Is unity." Thus, 
by definition, the characteristic largest value for a mission of length T 
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is determined from Equation (gti) by setting the expected number of threshold 
crossings to unity. The required characteristic largest value for the stand- 

I 

ardized normal variate is then 

G = 12 «n(T-E[H^(0)])]’^2 

where o is defined by Equation (59) and 

Elli^(0)3 is defined by either of Equations ( 6 l) or (62). 

This characteristic largest value for a stationary Gaussian random process 
having zero mean is sufficient to completely define any of the three dis- 
tributions used for normal extremes. 


As discussed in Section 2.1, the convergence of normal extremes to the type 
I asymptotic extreme-value distribution is extremely slow. The type 1 dis- 
tribution is therefore recommended for describing normal extremes only v/hen 
the characteristic largest value for the standardized variate exceeds 8. The 
extremal intensity function (q) corresponding to the standardized character- 
istic largest value (u) is given by Gumbel (Reference 1, page 137) as follows 
for normal extremes: 


a = u + u ^ - 2 u ^ + 


lOu 


-5 


(27) 


The extremal type I cumulative distribution function for the standardized 
variate (y) is given by 

«’^*^(y) = exp(-e “^) ( 26 ) 

where u and a are defined by Equations (65) -'nd { 27 ), respectively. 

The extremal type III distribution, v/hich converges to the type I distribu- 
tion v/ith increasing u, is recommended for characteristic largest values 
betv/een 3 and 8 for the standardized normal variate. The type III distribu- 
tion function for normal extremes was first suggested by Fisher and Tippett 
(Reference 5)- A sui table form of this distribution for the standardized 
variate is 

= exp|^-{ ^ j ( 23 ) 


RTTHODUCi::I'.lIV OF THE 
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where k = 




{u^-D 


(2A) 


The percentage points of this distribution as a function of the curaulative 
probability, p, are given by Gurbel {Reference 1, page 299) as 



.r. 


exp|«.n u 


2n(-£no) 


L 


] 


(25) 


As discussed In Section 2.1, for characteristic largest values of the 
standardized normal extremes close to 2, the lognornal distribution is 
essentially identical to the actual distribution of normal extremes cal- 
culated by Tippett and plotted in Reference 1 (page 129). For values less 
than 3» the lognornal approximation is generally more accurate than the 
extremal type III aoproximation and is therefore recommended for this ranee. 
The lognormal probability density function Is 


f(x) =— L-exp- i ( (16) 

^ «x 2 ° 

v/here y = tnx+ino (66) 

6= Un(J)]’''2 (67) 

X 

X is the median of the standardized normal extreme, and 
X is the mode of the standardized normal extreme. 

The required median (x) of the standardized normal extreme is obtained from 
the following equation: 


F('3^) = exp[0.63315(F(u)-l)] 


(20) 


where F(u) is the normal cumulative distribution function evaluated at u. 
Equation (20) is a modified fc-m of the equation for medians of extreme values 
given by Gumbel (Reference 1, page 79)- The Gaussian probability functions 
are tabulated, for example, in Reference 8 (page 555) and Reference 24 
(page 33) • 

The required mode (x) of the standardized normal extreme is obtained from the 
follo./ing equation: 
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x*F(x) F(u) (22) 

f(x) l-F(u) 

Equation (22) is a modified form of the equation for modes of normal extremes 
given by Gumbei (Reference ] , page 133)- Tabulated values of F(x)/f(x) 
may be found in Reference 2 (page 11). 


33 



NUMERICAL DEMONSTRATION 


ij.O 

The me*’hods previously described for determining limit-load probability 
distributions from time-domain and frequency-domain analyses have certain 
limitations which may be best illustrated by numerical exanples. Section 
^.1 presents three numerical examples of the Taylor's series method which 
demonstrate the effects of the method's fundamental assumptions. Section 
h.2 demonstrates the statistical estimation method which may be used to 
reduce the required number of Monte Carlo simulations. These numerical 
examples are based on sets of random numbers generated by a digital computer. 
The method for determining limit loads from a frequency-domain analysis is 
demonstrated using numerical data obtained from an analog Gaussian noise 
generator* Section k.3 presents the results of three examples of this 
method. 




I * 




Examples of Taylor's Series Method 


As discussed in Section 3*1 i the Taylor's series method for estimating the 
probability distribution of a nonlinear function of several random variables 
is based on the following three assumptions: 

(1) that the higher-order terms in the Taylor expansion- of the function 
are negligible compared with the first order terms, 

(2) that the Individual random variables are mutually independent, and 

(3) that the individual random variables are each normally distributed. 

The following is a brief discussion of the implications of these assumptions 
with numerical examples. 

Consider a function of four random variables 


f(W,X,Xi} = 1/V' X/Y + Z 

^ ^ = 7.J Sk O.S / / 


where 


Z"' -5o) 

By the Taylor's series method, the estimates of the mean, variance, and 
standard deviation of the function are as follows: 
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The stochastic behavior of this function was studied for three different 
cases. Case 1 involved dependent variables “ "0*5) af’d a non-normai 

variable with the variable X being uniformly distributed in the range 1 
to 3- Case 2 involved a non-normal variable (X n, U (1,3)) hut all variables 
were independent. Case 3 involved all normal and independent variables. 

The mean and standard deviation and the cumulative distribution function 
(CDF) were determined from a Monte Carlo simulation using a sample size of 
2000 for each of the three cases. The Monte Carlo simulations v;ere performed 
with the Boeing Generalized Statistics Program (GESP) described in References 
13 and 1^. The resulting means and standard deviations are presented In 
Table ^.1-1 for comparison with the Taylor's series estimates. Results of 
significance tests of the hypothesis that the Monte Carlo parameters are 
Identical to the Taylor's series parameters are also presented in Table 
l(.l-l along with the results of a Chi-square test for normality (Reference 
8, page 366). The hypothesis test for the mean was performed using Student's 
t statistic (Reference 8, page 127). The hypothesis test for the standard 
deviation was performed using the Chi-square statistic (Reference 8, page 138). 
The acceptance probabilities for such hypothesis tests are usually established 
at either c:.e percent or five percent levels. Values of the 37”degree-of- 
freed'-m Chi-square statistic corresponding to these probability levels are 
59 and 52, respectively. 

For this particular function, the Monte Carlo means and standard deviations 
are seen to approach the Taylor's series parameters as the assumptions of 
independence and normality of the individual variables are better satisfied. 

The hypothesis tests indicate that the mean determined by the Taylor's series 
method is sufficiently accurate regardless of normality and independence of 
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the individual random variables; the standard deviation determined by the 

I 

Taylor's series method is sufficiently accurate only when the individual 
random variables are independent. However, for none of the three cases 
was the hypothesis of normality verified by the Chi-square test. In Figure 
^.l-l, the cumulative distribution function determined from the Monte Carlo 
simulation for Case 3 is plotted versus the Taylor's series normal distribution 
to illustrate the results of the Chi-square test. 

These numerical examples are consistent with the theory discussed in 
Section 3.1. An accurate estimate of the mean requires only that the 
higher-order terms in the Taylor's series expansion are negligible, v/hereas 
an accurate estimate of the variance requires the additional assumption of 
independence among the individual random variables. All three assumptions 
must be satisfied in order that the function be approximately normally 
distributed. For the function studied, the second and higher partial 
derivatives are negligible or zero except with respect to the Y variable. 

The numerical influence of the neglected non-zero terms on the Taylor's 
series estimate of the mean and standard deviation appears to be small. 
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k.Z Examples of Monte Carlo Method 

As described in Section 3.2, the statistical estimation methe;i is a technique 
for conservatively estimating the lognormal limit-load parameters from a 
small sample of observed loads from a Monte Carlo simulation. The numerical 
demonstration of this method is based on a simulated analysis in which the 
limit load is defined as the largest load occurring in 100 independent 
observations of a standardized normal variate. By means of the GESP random 
number generator (References 13 and l4), 2000 simulated limit loads were 
generated. The limit-load distribution was approximately lognormal 
(P[x 37 > ^9l = 0.095) with mean equal to 2.509 and coefficient of variation 
equal to 0.1715. Ten data sets of 20 values each were statistically analyzed 
to determine the sample mean of the logarithms (m*) and the standard deviation 
of the logarithms (s*) . The best-fit mean m {i>) and rtandard deviation s (2) 
of the logarithms were then conservatively estimated using the 90 ^ one-sided 
confidence limit (Equation 51) for two sets of probability levels. The 
probability levels designated confidence fit ‘‘a” v/ere biased to positive values: 
= 1, 2, 3. 5. The probability levels designated confidence fit "b" 

were unbiased: = -k, ~3, -2, -1, 0, I, 2, 3i 5« The sample data and 

the conservative estimates for confidence fits *'a" and *'b" are presented in 
Tables A. 2-1 and 4.2-2, respectively. For comparison purposes, the “true" 
sample mean and standard dr-/iation of the logarithms based on 2000 values 
are = 0.505 and s^ = 0.1703. 

The data presented in Tables 4.2-1 and 4.2-2 are plotted on normal probability 
paper in Figures 4.2-1 through 4.2-10. Each plot shows, for each data set, 
the conservatively estimated distributions based on 20 values along with the 
“true" distribution based on 2000 values. Both conservative distributions 
result in values larger that the "true" values for the probability range of 
interest. Values from the biased confidence fit "a" suggest that most of 
the conservatism is in the estimate of the standard deviation. Values from 
the unbiased confidence fit "b" show a more balanced approximation to the 
"true" distribution. 
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Table k. 1-1 Parameters for Statistical Estimation 
Demonstration Using 90% Confidence Fit "b" 


Data Set 

1 m* 

- V 

s* 

y 

m (6) 

Sy (B) 

1 

1 0.886 

^ -11 

0. 1830 f 

1.004 

0.1995 

2 

0.80i» 

0.1669 : 

0.912 ! 

0.1820 

3 

0.893 

0.2320 T 

1.043 I 

0.2530 

k 

0.872 j 

0.1888 : 

0.994 1 

0.2058 

5 

0.910 j 

0.1319 I 

1.027 j 

0.1983 

6 

0.961 

0.1459 [ 

0.995 i 

0.1591 


0.877 T 

0.1467 i 

0.972 1 

0.1600 

8 1 



0 

0 

_ 0.1543 j 

1.030 i 

0.1682 

9 i 

0.815 i 

0.i684 . 

i 

0.924 ; 

0. i^~‘ 

10 1 

0.859 ■ 

0.1750 ^ 

0.971 
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^’3 Examples of Frequency - Domain Method 

Numerically demonstrating the method for determining tiie limit-Ioad 
probability distribution from the power spectral density function (PSD) 
of a Gaussian random process has two general aspects. The first is 
demonstrating that the extreme values from a continuous Gaussian time 
series of specified duration behave mathematically as extremes from a 
population of discrete normal variates. The second is demonstrating that 
the lognormal and extremal type 111 distributions provide valid repre- 
sentations of the actual distribution of extremes when based on the follow- 
ing expression for the standardized characteristic largest value derived 
in Section 3.3: 

[2 tn(T-E[N^(0)l)l*^^ (65) 

Both of these aspects will be demonstrated using numerical data obtained 
from an analog Gaussian noise generator. 

The Elgenco Model 311A Gaussian Noise Generator v/as used to obtain the 
required random time histories. This electronic device provides a stable 
and reliable source of Gaussian random noise having the following 
characteristics: 

(1) The output PSD is uniform to +0.1 dB from 0 to 35 Hz; the output 
falls off rapidly above ifO Hz. 

(2) The amplitude probability density function is Gaussian (normal) to 
less than +1 percent. 

The output of the Gaussian noise generator was passed through three first- 
order filters, all having cutoff frequencies of 25 Hz. Tine purpose of 
this filtering v;as to specify accurately the high-frequency roll-off so 
that the actual PSD could be precisely defined. The PSD used for the 
numerical demonstration is defined as 
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r(f) 


for 0 < f < ^0 


( 68 ) 


a^f S 

c 


2^3 




for f > i|0 


where f^ = 25 Hz, and 
2 

a = 2.785 <s the magnitude factor determined 
empirically from the generated output- 


Figure 4.3"! presents time histories of the unfiltered random noise pro- 
duced directly by the Elgenco Hoise Generator and of the random noise after 
it was passed through three 25 Hz filters. 

To determine the actual mean and standard deviation of the generated time 
series, a statistical analysis of the filtered Output time history was 
performed, based on the assumptions of ergodicity and stat ionari ty . A 
twenty-second duration of the output from the noise generator was sampled 
at 0.02-second intervals to provide 1000 data points. The mean and standard 
deviation of this large sample were then computed with the following 
results: 

V = -0.229 

o “ 6.418 

These statistical estimates were assumed to be the true parameters of the 
generated time series for all subsequent studies. 

A Chi-square goodness-of-f i t test (Reference 8, p. 3&5) was also performed 
with the sample of 1000 data points to verify that the generated output 
was Gaussian. The Chi-square statistic, based on a division of the data 
into 19 cells, was 19. This value corresponds to a Chi-square cumulative 
probability of less than 75 percent. Therefore, the random time histories 
obtained from the noise generator may be justifiably considered Gaussian 
with parametric values as estimated. 






The magnitude factor of 2.785 used in the PSD expression defined by 
Equation (68) is consistent -with the empirically determined standard 
^ deviation. With this PSD expression, the expected rate of zero crossings 

per second with positive slope va3s calculated to be 12.77. This compares 
well with the observed average of 13.0 taken from 26 seconds of the generated 
random time history. 

Because the theory developed in Section 3-3 applies only to Gaussian random 

processes having zero mean values, the DC bias of -0.229 was subtracted 

from the generated time histories for all subsequent numerical comparisons. . 

Each desired extreme-value observation was defined as the largest positive 

value occurring in a specified time duration, T. The theoretical symmetry 

of the Gaussian distribution was used to assure uniformity of the extreme- 

value data. Thus, observations of extremes were taken equally from the 

positive and negative peaks, and the absolute values of the tvxj data sets 

were combined into one total sample. The positive and negative data values 

were selected from different sections of the random time history, except 

for the data corresponding to T = 100 seconds. According to Gumbel 

(Reference 1, p. 110), the extreme largest and extreme smallest values are 

asymptotically independent for large samples. Thus, the combining of 

positive and negative extreme values from the same time-history section for 

T » 100 seconds is believed to be justified theoretically. 

Verifying that observed extreme values from a generated time series behave 'i 

mathemat ical ly as normal extremes v/as accomplished by comparing observed 

i 

and theoretical cumulative probability distributions. The theoretical | 

probability disti ibut ion, representing the largest Individual in samples j 

of size n taken from a standardized normal population, vas tabulated by ^ 

K. Pearson in Reference 2, (page 162). The basic relation betv/een sample | 

size (n) and characteristic largest value (u) for a specified distribution 
is given by Gumbel (Reference 1, page 82) by | 

F(u) = I - i (2) 

n 

I 
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The theoretical distribution of standardized normal extremes for 200 dii- 
crete samples, which corresponds to a standardized characteristic largest 
value of 2 . 5758 , was selected for comparison wi th an observed extreme- value 
distribution obtained from the generated time series. From Equation (65), 
the value of T corresponding to this character istic largest value is 2.161 
seconds for the specified PSD. The observed probability distribution v;as 
obtained from 200 samples of largest values occurring in time intervals of 
2.161 seconds from the generated random time history. The 200 values were 
modified to eliminate the DC bias and then ranked in increa'*’'^ order. The 
cumulative probability assigned to the ith observed value v.'as i/201. 

The comparison betv;een the observed and theoretical disti '.'ions shown 
in Figure f|.3-2. The excellent agreement is corroborated lj . ( s-jare 
test of the hypothesis that the observed distribution is identical to the 
theoretical distribution of normal extremes. The Chi-square statistic, 
based on a division of the data into 20 cells, was 22.8. This value 
corresponds to a Chi-square cumulative probability of approx iinately 80 
percent. Therefore, the observed extreme values obtained from the generated 
random time series may be considered as normal extremes from a discrete 
sample of size n, v;here n is determined from Equation { 2 ) given the 
characteristic largest value. 

Verifying the accuracy of the lognormal and extremal type 111 distributions 
to represent normal extremes in terns of the standardized characteristic 
largest value from Equation ( 65 ) v^as accomplished by comparing the approxi- 
mating distributions vnth the distribution of observed extremes for tvK> 
different time intervals, T, For T = 1.0 second, the standardized 
characteristic largest value corresponding to the specified PSD is 2.257. 

This is in the range where the lognormal distribution provides a nearly 
perfect representation of the actual distribution of normal extremes. 
Therefore, the required lognormal parameters may be obtained from Equations 
(20), (22), ( 66 ), and (67) of Section 3 . 3 . The Straight line in Figure l*.3-3 
corresponds to y = 2.733 and 5 = 0.193- The lognormal approximation is seen 
to provide a very good representation of the observed distribution plotted 
from 500 data points. 
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For T = 100 seconds, the standardized characteristic largest value 
corresponding to the specified PSD is 3-782. This is in the range v/here 
the extrenal type III distribution provides a somewhat better representation 
of the actual distribution of normal extremes than the lognormal distri- 
bution. The required extremal type III parameter for the standardized 
variate, obtained from Equation (2^*) of Section 3.3, is k = 17-60A. The 
corresponding distribution, obtained fro.m Equation (25) of Section 3-3, is 
shown in Figure The lognormal distribution also shown in this 

figure corresponds to y = 3-213 and 6 = 0. 1*»0. The extremal type III is 
seen to provide an adequate representation of the observed distribution plotted 
from 118 data -points. The standard deviation for the lognormal approximation, 
which determines the slope of the straight line, is seen to be too large. 

Sinca the median value is accurately determined for both approximating 
distributions, the character istic largest value calculated fiom the Gaussian 
PSD by Equation (65) is seen to be the proper value for defining the desired 
extreme-value distribution. 
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5.0 


CONCLUSIONS AND RECOKENDATIONS 


The methods for determining limit-load probability distributions from time- 
domain and f requency-doma in dynamic loads analyses have been described and 
numerically demonstrated. The primary contribution is obtaining the extreme- 
value probability distributions from the Gaussian PSD of a frequency-domain 
analysis- Another contribution is obtaining conservat iv» estimates nf the 
limit-load probability distributions from a small number of Monte Carlo 
simulations. 

Recommended areas for additional research include the following: 

1. Improve the accuracy of the estimate of the parameter 5 in Equation (67) 
for the lognormal approximation to the distribution of normal extremes 
v/hen u > 3- 

2. Study the possibility of combining the Monte Carlo technique of Russian 
Roulette and the statistical estimation method to determine more efficient 
estimates of the limit-load parameters from Monte Carlo analyses. 

3. Develop a method for obtaining the extreme-value distributions frcn a 
combination of time-domain and frequency-domain dynamic loads analyses. 

k. Study methods of accounting for payload mass and stiffness variations 
in dynamic loads analyses. 
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